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Abstract 

The chiral hmit k ~ ndP) in lattice gauge theories with Wilson fermions and 
problems related to near-to-zero ('exceptional') eigenvalues of the fermionic matrix 
are studied. For this purpose we employ compact lattice QED in the confinement 
phase. A new estimator for the calculation of the pseudoscalar mass m-Tr 

is proposed which does not suffer from 'divergent' contributions at k ~ Kc{i3)- 
We conclude that the main contribution to the pion mass comes from larger modes, 
and 'exceptional' eigenvalues play no physical role. The behaviour of the subtracted 
chiral condensate {ipip)suht near Kc{I3) is determined. We observe a comparatively 
large value of {ipip)subt ■ Zp^ , which could be interpreted as a possible effect of the 
quenched approximation. 

1 Introduction 

As it is well known, chiral symmetry is broken explicitly in lattice gauge theories 
with Wilson fermions as in QCD and QED [|l|. Presumably, it can be restored by 
fine-tuning the parameters in the continuum limit. If so, one can approach the 
continuum limit and chiral symmetry restoration along a 'critical' line /tc(/3). It 
is another question, whether on this line the chiral symmetry becomes explicitly 
realized or spontaneously broken. In the continuum limit the lowest-lying state in 
the spectrum of the Wilson fermionic matrix M.{k,; U) should have eigenvalue zero 

for K — ^ KciP) ■ 

For nonzero lattice spacing a 7^ the chiral symmetry cannot be restored 
exactly. At the same time in the confinement phase the pion mass tends to 
zero in the limit k ^dP) (the so-called partial symmetry restoration). The 
mechanism of this partial symmetry restoration is still not well determined. For 
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example, one can not exclude that it is connected with the transition to a parity 
violating phase where (tpj^tp) 7^ [@]. If this transition is of second order then 
there will be a massless pion at k = ndP)- 

Computations near the chiral transition line Kc(/3) are notoriously difficult. 
The first problem is that the matrix inversion becomes very slow when k, ap- 
proaches Kc{P) and most of the known inversion methods fail very close to Kc(/5). 
Nevertheless, the conjugate gradient method Q appears to be reliable even in 
the 'critical' region. Another problem which arises in simulations of lattice gauge 
theories with Wilson fermions is that near the chiral transition very small ('excep- 
tional') eigenvalues {Aj ~ 0} of the fermionic matrix A^(/t;f/) appear which 
make it practically impossible to approach ndP) without an enormous increase 
of statistics. The usual way to bypass this problem is to carry out calculations at 
K-values sufficiently below ndP) and then to extrapolate the observables, e.g., 
hadron masses, to the 'critical' value Kc{(3). 

Apart from practical considerations the zero mode problem is connected with 
the question of the mechanism of the chiral transition in lattice theories with Wilson 
fermions. As it was already mentioned, zero eigenvalues are expected to exist in the 
chiral limit in the continuum. Within this context the question has to be answered, 
whether the near-to-zero eigenvalues of A4{k,; U) observed on the lattice in the 
'critical' region in finite volumes have physical relevance or rather arise as an effect 
of the lattice discretization. 

Configurations with extremely small eigenvalues {Xi} were first discovered 
in P, H, and were called 'exceptional' configurations. In fact, very close to the 
transition point k ^ ndP) such configurations appear to be very 'normal' (see, 
e.g. []^, ^ 1^, )• In practice the appearance of these configurations can be used as 
an indicator for approaching ndP)- 

It is important to point out that the appearance of near-to-zero eigenvalues 
{Xi} at K, ^ ndP) is not a disease of the quenched approximation. The fermionic 
determinant decreases the spread of small eigenvalues but does not eliminate them 
totally i,0,|,|ig. 

It is the aim of this work to study the properties of a lattice gauge theory with 
Wilson fermions very close to the chiral transition line ndP). In the confinement 
phase we study the behaviour of the pseudoscalar pion mass uIt, and the subtracted 
chiral condensate {ipip)subt , which for Wilson fermions is an order parameter of 
chiral symmetry breaking in the continuum limit, near ndP). We employ compact 
lattice QED in the quenched approximation. Compact lattice QED with Wilson 
fermions in the confinement phase possesses similar features near the chiral tran- 
sition K K,dP) nonabelian theories, and thus can serve as a comparatively 
simple model of QCD in this study. The cost of numerical calculations is much less 
than for QCD. We expect that our analysis is applicable to QCD, as well. 
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2 Action and observables 

The standard Wilson lattice action Swa{U, ip, ip) for id compact U{1) gauge theory 
(QED) is 

Swa = Sg{U)+Sf{U,^,^). (1) 
Ineq.(P Sg{U) is the plaquette (Wilson) action for the pure gauge U{1) theory 

SGiu)=p-J2ii-cosep) , (2) 

p 

where /? = l/g^are ' UxfM = exp{i9x^), 6^^ G (— 7r,7r] are the field variables 
defined on the links / = . Plaquette angles 6p = O^.^^, are given by O^.^^i, = 

The fermionic part of the action SpiU^ip^ip) is 



x,y s,s' = l 



M 



i-K-Q{U), 



Q 



xy 



E[^: 



y,x+fi 



(3) 



where Ai is Wilson's fermionic matrix, and k is the hopping parameter. 

It is known, that at least in the strong coupling region for the standard Wilson 
action (in the confinement phase) the ordinary mass term and the Wilson mass term 
cancel within the pseudoscalar mass m^r at some k = ndP) , so that quadratic 
terms in the effective potential vanish for the pseudoscalar field , and Hc ~ 0.25 at 



/? = 12 



In the weak coupling range perturbative calculations indicate that the mass of 
the fermion becomes equal to zero along the line Kc{f3) (for the free field theory 
= 0.125) |TT[. 

We calculated the following fermionic observables 



(n) = -L.(Tr(A^-S5-M-S5))G, (4) 

where ( )g stands for averaging over gauge field configurations, and V = Nr- 
is the number of sites. Pseudoscalar zero-momentum correlators r(r) are defined 
as follows 
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r(r) = --^ ■ ^ (V'75^(^,^) ■ V'75V'(0,y)) 

s x,y 

where Sp means the trace with respect to the Dirac indices. These correlators 
as well as the pion norm 11 appear to be very sensitive observables in the 'critical' 
region. This can be understood by considering the spectral representation of the 
fermionic order parameters. Let /„ = fn{s,x) be the eigenvectors of Ai with 
eigenvalues A„ , and gn = gn{s,x) be the eigenvectors of ■y^M. with eigenvalues 

Mfn = \n- fn , I^MQu = fJ'n ' g-n. ■ (6) 

Then one can easily obtain a spectral representation of the fermionic order para- 
meters : 



Evidently, an eigenstate of with eigenvalue zero is also an eigenstate of ■j^M. . 
So, the presence of configurations which belong to zero eigenvalues of A4 also gives 
rise to poles in 11 . For correlators one obtains 

Y: Sp(A1;oS5-Moi75) 1.4=. = E — ■ bnnir) ■ — ■ bn'niO) , (8) 
X nn' f^^' 

where 

bnn'ir) = J291{s,x,t) ■ gn'{s,x,T) . (9) 



For further discussion on properties of the fermionic matrix see, e.g., [p!3| , ^ 
and references therein. 



3 The pion mass in the chiral Hmit 
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3.1 Pseudoscalar correlators and the standard estimator 
ttItt near Kc{f^) 

Transfer matrix arguments suggest the following form of the pseudoscalar correlator 
r(r) (at least, for t > ) 



T(t) = 



+ higher-energy states . (10) 

The standard choice of the estimator for the effective mass of the pseudoscalar 
particle m^-^(r) = m-^ij) is 



cosh m^(r) -r- l) _ f(r + 1) 
coshm,(r)(^-r) ~ f(r) 

= e-'^-(") if m^Nr > 1 , 
where r(r) is the estimator of the pseudoscalar zero-momentum correlator 

r(r) = --Er.(r) , (12) 

and n is the number of measurements. The r-dependence of the effective mass 
stems from the contribution of higher-energy eigenstates of the transfer matrix. 
With increasing r the contribution of these higher-energy states is expected to 
be suppressed, and the resulting plateau in the r-dependence of the effective mass 
gives the true mass itIt, . 

However this approach fails when one comes close to the chiral limit, i.e., when 
K, KciP) ■ The well-known problem in both QED and QCD with the calculation of 
(and other fermionic observables) is connected with extremely small eigenvalues 
of the fermionic matrix (and correspondingly 75 A1) which appear on a finite 
lattice at some k'{(3]Ns]Nt-)'^Kc{(3) . 

To illustrate the problem we show in Figs.la-c the dependence of the observables 
{'iljip) , {'ip'j^ip) and (n) on K at P = 0.8 on a 8'^ lattice. Well below the 
transition point, i.e., when k < k'^{[3] N^] Ng) , these averages are statistically 
well-defined. The increase of the number of measurements from, say, n = 100 to 
n = 200 produces just a shght change of the averages, and statistical errors decrease 
as ~ ■ The situation however changes significantly at k^k'^{(3] N^] Ng) (to 

the right of the vertical dashed line in Figs.la-c). The averages begin to fluctuate 
drastically with increasing n (compare circles, crosses and squares in Figs.la- 
c), and the errorbars become dramatically large. It is worthwhile to note that 
the increase of statistics does not necessarily entail the diminishing of the errorbars 
(compare, e.g., circles and squares in Fig.lc). This means in fact that these averages 
and errors, which were calculated using the jackknife procedure, make not very much 
sense and the accumulation of measurements will not change the state of affairs. 
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Figure 1: The dependence on k of {i'i') (a), {^Jl^^j) (b) and (IT) (c) 
(3 = 0.8 on a 8^ lattice. Dotted lines have been added to guide the eye. 



at 



Note, that Nr] Ng) in Fig.lc is a little bit smaller than in Figs.la,b, since 

(n) is the more sensitive observable with respect to small eigenvalues of 75AI 
(and M. respectively). 

The behaviour of (ip'y^ip) in Fig. lb deserves maybe some additional comment. 
It was suggested [0] that at k, = ndP) there is a transition to a parity- violating 
phase, in which {4''j5ip) 7^ (see also [|r§]). At first sight our data for (V'TsV') 
in Fig. lb look like the confirmation of this hypothesis, but in fact they are not, 
because all averages and statistical errors to the right of the vertical dashed line are 
unreliable, and no definite conclusion can be drawn. 

It is interesting to mention that on 'exceptional' configurations the number of 
conjugate gradient iterations Ncg behaves differently from the 'non-exceptional' 
case where the convergence of the CG method is mainly determined by the condition 
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Figure 2: The time history of 11 and Nf.g at j3 
S'^-lattice. 



0.8 and k = 0.22 on a 



number ^ = Xmax/^min of the to be inverted matrix [|. Moreover, when Amm is 
sufficiently small then Ncg behaves as 



_ 1 

Ncg ~ ^min 



(13) 



and therefore can be used as an indicator for approaching the chiral limit (see, e.g., 
P, Q). However, it is generally not known for which classes of distributions of 
eigenvalues eg. ([T3|) holds. In a special study of this problem using hermitian ma- 
trices with some predefined (but continuous) distributions of eigenvalues we could 
confirm the commonly accepted results |jlO|- Now, in case of 'exceptional' small 
eigenvalues the behaviour of Ncg does not follow eg. (pTS]) . This can be seen from 
Fig.^], where we display the time histories of the pion norm 11 and Ncg on a 
8^ lattice slightly above Kc{I3). While the pion norm 11 develops huge spikes 



here A, 



max / vain 



denotes the maximal/minimal eigenvalue of the matrix to be inverted 
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(amplitudes up to ~ 10^) for some configurations, the values of N^g vary within 
a ~ 10% corridor around the average value. Thus the CG convergence behaviour 
seems to be not only determined by ^ but may also considerably depend on the 
eigenvalue distribution of the matrix to be inverted. 
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Figure 3: The r-dependence of the standard estimator 
lattice size is 12 ■ 4^. Lines are to guide the eye. 



nir. 



T 

(r) at (3 = 0; the 



In Fig. 3 we present the r-dependence of the standard estimator mT,{T) eq.(|TT]) 
for different k's at = on a 12-4^ lattice. The value of Kc is expected to 
be ~ 1/4 at this value of (3 (at least in the infinite volume). For the 12 ■ 4^ 
lattice the 'critical' region starts at k'^{(3 = 0) ~ 0.24. With increasing the lattice 
size this 'critical' region shrinks : Nr] Ng) Hc{P) for A^^, A^s ^ oo 0. 

Nevertheless, on a finite lattice the problem connected with the approach to the 
chiral limit remains, and the usual 'safe' choice of k in QCD calculations, i.e., the 
choice of k well below Kc(/9), entails an unrealistically large ratio niT^/mp . 
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Sufficiently below the 'critical' value of k the standard estimator of the pion 
mass eq. ([TlD yields well-defined results. The contributions of higher excited states 
die out with increasing r, leaving a nice plateau in the r-dependence of the effective 
mass m7r(T), which determines the actual mass at these values of couplings. For 
the reasons given above, this estimator fails to give precise results in the region 
^ > f^'ciP] N-r] Ns) , i.e. when approaching the chiral limit. (All data points of this 
figure represent averages of at least 5000 measurements). 



3.2 Another estimator for m 



IT 



One possible way to bypass this problem is based on the observation which we call 
the factorization of the contribution of near-to-zero eigenvalues. In Figs.4a-b we 
display the values of rj(r) on individual configurations i = 1, . . . , 900 on a 16-8^ 
lattice at /5 = and k = 0.25 for r = 3 (Fig.4a) and r = 4 (Fig.4b). The huge 
spikes (up to ~ 10'^) in the time histories are due to the existence of very small 
eigenvalues in the spectrum of •y^M. (respectively M.). These spikes appear with 
a certain frequency and just reflect the content of the spectrum of 'j^M. for this 
value of K, which is close to ^dP) in a finite system, i.e. the lowest eigenvalues are 
allowed to fluctuate from configuration to configuration. 

The important observation is the absence of such peaks in the ratio of the 
individual correlators 

«.(^) - (14) 



which means the factorization of the 'divergent' contributions originating from near- 
to-zero eigenvalues. The statistical fluctuations of this ratio are very small even at 
(/5) (see Fig. 4c) as compared with fluctuations of rj(r) (Figs.4a,b). 
Therefore, we propose to use the ratio eq.(|l4D for the extraction of the 'pion' 
mass near the chiral transition. Then a new estimator for the pseudoscalar mass, 
which we denote by m7r(T) , can be obtained from the following expression 

cosh m.(r) -r-l) ^ i ^ ^^^^ ^^^^ 
cosh m^(r) - n .^^ 

Therefore, the relation to the standard estimator m^(r) in eq.(p!lD is given by 



cosh m^(r) -r- l) cosh m^(r) - r - l) i sTAt) , , 

\. 7V(^r ^ = , - , .(N^ — ^ ■ ^ (16) 

cosh m7r(T) (-2^ — rl coshm,r(T)' — 



i=l 



Ffr) 



with 



5r,(r) = r,(r) - f(r) ; 5^5r,(r) = 0. 
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Figure 4: Time histories {rj(r)} for r — 3 (a), r — A (b) and 
{ri(T = 4)/ri(T = 3)} (c)ona 16-8^ lattice at /3 = and k = 0.25. 

In general, of course, this is not the correct way to calculate masses. However, 
under certain circumstances it is justified to do so. 

Assume that x and y are two correlated random variables with some dis- 
tribution P{x;y), where Jdxdy P{x;y) — 1. Then the average of any functional 
0{x;y) is 



(O) = Jdxdy O ■ P{x;y) . 
One can define the conditional average y{x) as 

y{x) = jdyyP{x;y)/ jdy P{x;y) . 
If the distribution P{x; y) is such that 



(17) 



(18) 
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r,(T=i) r,(T=i) 

Figure 5: rj(r = 2) as a function of rj(r = 1) at /3 = and k = 0.25 on a 
16 ■ 8^ lattice at different scales. Single points represent individual configurations. 
Straight lines correspond to the linear fit according to the least squares method. 
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(20) 



y{x) = C ■ X , C = const , (19) 

then it is easy to see that 

M = /I' 

(x) \x/ 

In our case x = rj(r) and y = Ti{T + 1). In principle, one can use also the 
ratios rj(r + A;)/rj(r) with k >2, though the signal-to-noise ratio decreases with 
increasing k. 

In Fig.|^ we show rj(2) as a function of at (3 = and k = 0.25 

on a 16 ■ 8^ lattice. For illustrative purposes we have shown our data of ~ 2000 
measurements with four different choices of the scale. The distribution gives a clear 
indication in favour of eg . (p!9D . The correlation coefficient p defined in a standard 
way il^ is very close to unity, and varies between ~ 0.98 and ~ 0.995 for 



different subsets of measurements shown in Fig.|^a -v- Fig.^. 

To provide a proof that fhT^{T) can serve as a reliable estimator of the pion 
mass mTj{T) we will compare the properties of both estimators. We'll show that 
for K-values sufficiently below Kc(/3) , i.e., where the standard estimator of m^r 
can be reliably defined, both estimators are in a very good agreement. For values 
of K very close to ndP), where the standard estimator fails to work, while 
fits the same straight line. 

We have run simulations in the confinement phase at (3 = and (3 = 0.8 
using lattices with A^^ = 12; 16; 20 and A^s = 4; 8 . For the matrix inversion we 
used the standard conjugate gradient method |] with even-odd decomposition 
TS| , |T^ which guaranties convergence even at k ^ i^c{(3)- 



The r-dependence of the effective mass m^(r) appears to be more complicated 
than that of the standard estimator m^(r) . In Fig. 6 we have plotted the depen- 
dence of ^^^(t) on T at different /?'s and k's on lattices with Ng = A and 
Nr = 12; 16; 20 . At (3 = the value k = 0.23 is sufficiently below the transition 
point /€c(0) ~ 0.25 (no appearance of 'exceptional' small eigenvalues of and 
75 respectively) and therefore ^^^(t) can be compared with the corresponding 
m7r(r) shown in Fig.3. The standard estimator m^(r) exhibits a nice plateau up 
to r = 6 on a 12-4^ lattice, while the corresponding plateau for m7r(r) is to 
be seen only for r^5 . In this case both estimators (m^(r) and m^(T) ) are in 
a very good agreement as long as r is not too close to A^r/2. These deviations 
from the standard m^(r) at r ~ A^r/2 are a general feature of the new estimator 
which presumably stem from the way of averaging the corresponding observables 
on a finite lattice (but which fade with increasing A^^, see below). The increasing 
of Nr entails the extension of the plateau. The same effect, i.e. the enlarging 
of the plateau in the r-dependence of rh.„{T) occurs for other /9's and k's, 
including values {(3) where the standard estimator is not well-defined (e.g., 

K, = 0.25; /3 = in Fig.6 ). It is important to notice that increasing A^^ does not 
change the position of the plateau, influencing only its extension. 
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Figure 6: The r-dependence of m7r(T) at = and (3 — 0.8 . Lattice sizes 
are 12 • 4^ , 16 • 4^ and 20 • 4^. Lines are to guide the eye. 



Finite volume effects, i.e., the dependence of rhniT) on Ng tend to increase 
shghtly with increasing couphng /3 - this could be interpreted as an extension of 
the 'critical' zone in k, with rising f3 . In Fig. 7a we compare the r-dependence 
of both estimators itIt^^t) and ttIt^^t) at f3 = and f3 = 0.8 for two lattice 
sizes : 16 ■ 4"^ and 16 ■ 8^ . The values of the corresponding /t's are chosen to be 
sufficiently far from the transition point, ensuring the applicability of the standard 
definition of the pion mass. At P — both estimators are in an excellent agreement 
even on the smaller lattice with A^^ = 4 . At P — 0.8 m^(T) on the smaller 
lattice gives an ?«3% overestimated value. However, on the lattice with Ns = 8 
the agreement between both estimators becomes very good. Note, that at these 
ft;-values the deviations of 171.^ (t) from the plateau for r — >• Nt/2 as described in 
Fig. 6 decreased substantially. 

In Fig. 7b and Fig. 7c we compare the values of rfiT^{T) at the same (3 's as in 
Fig. 7a but now at k's chosen near Kc{(3) ■ The finite volume dependence becomes 
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Figure 7: The r-dependence of m^(r) and rh-n-{T) at /9 = and /9 = 0.8 for 
two lattices : 16-4^ and 16 ■ 8^ (a); the r-dependence of ttIt^^t) at /? = 0.8 
(b) and /3 = (c) for lattices 20 • 4^ and 16-8^. Lines are to guide the eye. 

stronger when k becomes close to Kc{P)- The estimate of m^(r) is lowered by 
increasing the spatial extension Ng . The quahty of the plateau is good enough to 
determine m7r(T) with good accuracy. 

In Fig. 8 we show the ^-dependence of rfi^ for /3 = and P = 0.8 on a 16 ■ 8^ 
lattice. Different symbols correspond to different values of r. At both values of j3 
the new estimator behaves in a similar way. Even at >^c{l3) the values of 
at different r and /? are in a good agreement. For both /3- values shows 
a linear dependence on k in the vicinity of Kc{P). This leads to the behaviour 

ml = BM -m^-. m, ^ , (21) 
where the dimensionless bare fermion mass ma is defined as 



14 



, 6 



0.4- 



. 2 




19 .20 .21 .22 .23 .24 25 .26 



K 



Figure 8: The K-dependence of ml.{n) at /? = and P = 0.8 on a 16-8^ 
lattice for several values of r. The errorbars are smaller than the symbol size. 
Lines drawn are to guide the eye. 
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m„ = 7-^ . (22) 

At /3 = the straight line extrapolation (broken line in Fig. 8) predicts for the 
'critical' value Kc(0) ^ 0.2502(1) which is in a very good agreement with the strong 
coupling results |Tl|. With increasing /? the value ndP) decreases : Kc(/3 = 
0.8) ^ 0.2171(1) (errorbars correspond to one standard deviation in x^/'^d.o.f. )■ 
The corresponding slopes are : Bt,{0) = 4.91(4) and 5,^(0.8) = 3.42(3). 

For K-values sufficiently below ndP) , i-e., for n < tt'dP) , ml is just the same 
as the standard estimator . For k's very close to Kc{/3) the standard estimator 
can't be applied, still fits very well the same straight line having a very small 
statistical error. Our data at /3 = suggest, that for k ^ 0.246 (m^ ^ 0.034) the 
standard estimator of deviates from the straight-line behaviour in Fig. 8. For 
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instance, at k > 0.248 {rriq ^ 0.0176) the standard estimator was not applicable 
any more. Numbers change when considering the case /3 = 0.8. Here the deviation 
of from the straight-line behaviour starts already at rrig ~ 0.078, which 

corresponds to k = 0.21 . This demonstrates the advantage of 7?!^^ which gives 
the possibility to approach the chiral limit much closer than the standard estimator 
of m^r would allow. Especially, systematic errors induced by the finite volume can 
be investigated in the 'critical' region by means of rfi-jr with good accuracy, which 
in turn should allow a more precise extrapolation to the thermodynamic limit. 



4 The subtracted chiral condensate 

One can use the Goldstone theorem to define the 'physical' (subtracted) chiral 
condensate {tl'tl') suM ■ Following pO[ we define {ipip)subt which is supposed to 



serve as an order parameter for the chiral symmetry breaking (at least, in the 
continuum limit) 



— 11 — — 

(-ip-ip) subM = ( ) ■ ZpY.ii'l^'^^x ■'ip'^M , (23) 

where Zp is the renormalization constant of the pseudoscalar quark density and 
the sum in the r.h.s. of eq.(p3|) is connected with the pion norm (11) . At nonzero 
lattice spacing the pion mode in the limit /t — > Hc{l3) should result in a constant 
contribution to the condensate {'ip'ip)subt while the contribution of the higher energy 
(massive) modes should decrease as ~ (/tc — /t) . 

Fig. 9 shows the dependence of {'ijjil))subtiK-)- Zp^ on k for /5 = and 13 = 0.8. 
Here we have chosen k values in the region where near-to-zero 'exceptional' modes 
not yet show up. The pion mass reached already the asymptotic regime 

shown in Fig.8. For these k's the values of the condensate are on the (expected) 
straight line which permits to make a reasonable extrapolation to k — Kc{f3) . The 
extrapolated values of {ipip) subt{i^) are nonzero, of course (at least, because of the 
finite spacing). 

The main observation here is that the values of the condensate are comparatively 
large, comparing with what one could expect from the naive definition {ipil')bare — 
{ipip)pert with (ipip) bare ^ (although oue cauuot exclude that the reuormalizatiou 
constant Zp is responsible for this effect). 

A possible interpretation of this effect of amplification of {ip'4')subt ■ Zp^ is 
due to the unphysical contributions in the quenched approximation. In the dy- 
namical fermion case the fermionic determinant tends to decrease the pion norm 
substantially comparing with that in the quenched approximation JlOl (at least in 



the confinement phase). Therefore, we expect that in the dynamical fermion case 
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Figure 9: The mg-dependence of the subtracted condensate Zp • {i'i') subt{f^) at 
/3 = and P = 0.8 on a 16 ■ 8^-lattice. Lines are to guide the eye. 

the condensate should be much smaller. From another point of view the unphysical 
contributions in the quenched approximation were discussed also in ^ p^] . 

At K ~ Kc(/5) the behaviour of {i^i^) subtif^) should be mainly controlled by 
the value of the preexponential factor in the pseudoscalar correlator r(r) in 
eg. ([101). The existence of the pion requires to behave as ~ ^/ttl.^ when 
m^r — ^ . Together with eq.(pTl) this implies 

A^ ~ , m, ^ 0. (24) 

In Fig. 10 we present the value ' Tflq clS 3i function of ruq obtained from our 
data on a 16 ■ 8^ lattice. For bare masses niq < 0.06 we are not able to indicate 
A^ ■ rriq, since the extracted A^^'s become statistically unreliable. At (3 = the 
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value A^. ■ rriq shows a gradual increase with decreasing rriq, which, however, 
not in contradiction to eq.(^), taking into account the comparatively large values 
of TTiq. Remarkably, we observe A^. ■ rriq to fit a constant value at /3 = 0.8 with 
good accuracy in agreement with eq.fEil 
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at p = and p = 0.8 on a 16 ■ 8^ lattice. Added 



5 Conclusions 

In this work we have studied the chiral limit of a lattice gauge theory with Wilson 
fermions within the quenched approximation, employing compact lattice QED in 
the confinement phase. 

We observed that ratios of the pseudoscalar correlators Ti^r + l)/ri(r) do 
not suffer from near-to-zero ('exceptional') eigenmodes of the fermionic matrix. 
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This means that for every given configuration all 'divergent' contributions to the 
correlators are factorized. 

Making use of this observation we propose another estimator of the pseu- 
doscalar mass, which is well defined near the chiral transition line Kc{f3) and even 
for K > KciP) in contrast to the standard estimator mj^. 

It should be stressed that the possibility to introduce this estimator is based on 
the special kind of correlations between Ti^r) with different r , i.e., such that 
eq.(|l9D is fulfilled. This is not a universal property, of course. 

For K-values sufficiently below K-dP) , i.e., where the standard estimator of m^r 
can be reliably defined, both estimators are proved to be in a very good agreement. 
By approaching Kc{i3) we observe a linear dependence of on k : 



For values of k very close to ndP) the standard estimator fails to work, while 
still fits the same straight line having a very small statistical error. 

This study leads us to the conclusion that the main contribution to the pion 
mass comes from larger modes, and the 'exceptional' near-to-zero eigenvalues of the 
fermionic matrix 'y^Ai (and Ai respectively) play no physical role. This conclusion 
is in agreement with | p^4| . 

The new estimator gives the possibility to approach the chiral limit much 
closer than the standard estimator of would allow. Especially, systematic errors 
induced by the finite volume can be investigated in the 'critical' region by means of 
rh-jr with good accuracy, which in turn should allow a more precise extrapolation 
to the thermodynamic limit. 

The 'critical' value can be determined with high accuracy. Our ^^(O) 

is in very good agreement with strong coupling predictions 

We investigated the subtracted chiral condensate {4'4')subt following the 
definition in near the chiral transition line Kc(/3). We obtained a comparatively 



large value of {ipilj)subt ■ Zp which we interpret as a possible effect of the quenched 
approximation. 

It would be very interesting to know in which way the fermionic sea -dynamical 
fermions- could infiuence the behaviour of the pseudoscalar mass m-Tr and the 
condensate {ipip)subt ■ This work is in progress. 
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